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ABSTRACT 

We have conducted VLB A phase- referencing monitoring of H2O masers around the red su- 
pergiant, S Persei, for six years. We have fitted maser motions to a simple expanding-shell model 
with a common annual parallax and stellar proper motion, and obtained the annual parallax as 
0.413±0.017mas, and the stellar proper motion as (— 0.49±0.23 mas yr""'^, — 1.19±0.20 mas yr^^) 
in right ascension and declination, respectively. The obtained annual parallax corresponds to the 
trigonometric distance of 2.42l'^QQg kpc. Assuming the Galactocentric distance of the Sun of 
8.5 kpc, the circular rotational velocity of the LSR at the distance of the Sun of 220 km s~^, 
and a flat Galactic rotation curve, S Persei is suggested to have a non-circular motion deviat- 
ing from the Galactic circular rotation for 15 km s^^, which is mainly dominated by the anti 
rotation direction component of 12.9 ± 2.9 km s~^. This red supergiant is thought to belong to 
the OB association. Per OBI, so that this non-circular motion is representative of a motion of 
the OB association in the Milky Way. This non-circular motion is somewhat larger than that 
explained by the standard density-wave theory for a spiral galaxy, and is attributed to either 
a cluster shuffling of the OB association, or to non-linear interactions between non-stationary 
spiral arms and multi-phase interstellar media. The latter comes from a new view of a spiral arm 
formation in the Milky Way suggested by recent large N-body/smoothed particle hydrodynamics 
numerical simulations. 

Subject headings: stars: supergiants, masers, Galaxy: structure 



1. Introduction 

Phase-referencing VLBI technique is capable of 
determining positions of celestial objects relative 
to background calibrators with sub-milliarcsecond 
accuracy (Bartel et al. 1986). Absolute position 
changes of Galactic maser sources can be directly 
detected on the sky using this technique, and 
the observables are very useful for measuring an- 
nual parallaxes or trigonometric distances to the 
sources and their proper motions (e.g., Xu et al. 
2006; Hachisuka et al. 2006). In addition, the mea- 
sured trigonometric distance and proper motion 
together with the radial velocity of a source can 
provide us a complete three-dimensional kinematic 
information without model assumptions. 

There are a number of Galactic evolved stars 
with circumstellar H2O masers, which are observ- 
able with VLBI at 22 GHz. The red supergiant, 
S Persei (hereafter, S Per), is one such target har- 
boring H2O masers as well as SiO and OH masers 
(e.g., Alcolea et al. 1990; Richards 1996). This star 
is separated by 1.5° from a double cluster, h and 
X Persei (hereafter, h + x Per, see Figure 1), lo- 
cated at the outer edge of the Perseus arm, which 
involves a large group of more than 20 red super- 
giant stars. The cluster is the most massive open 
cluster within 3 kpc from the Sun and belongs to 
the OB association. Per OBI. Humphreys (1975) 
reported that the radial velocity of S Per is fairly 



consistent with that oi h + x Per, so that S Per 
is thought to belong to the OB association. We 
can therefore determine the distance and three- 
dimensional motion of the OB association in the 
spiral arm by measuring the annual parallax and 
proper motion of the member star. 

We report results of phase-referencing VLBI 
monitoring of the H2O masers associated with 
S Per over six years with the Very Long Base- 
line Array (VLBA) of the National Radio Astron- 
omy Observatory (NRAO). The observations are 
described in Section 2. Data reduction including 
astrometric analyses is presented in Section 3. The 
results are presented in Section 4. We discuss the 
results in Section 5 and summarize this study in 
Section 6. The results presented here replace those 
reported previously by Asaki et al. (2007a). In 
this paper, we define that the Local Standard of 
Rest (LSR) is a reference frame at the position of 
the Sun and moving in a circle around the cen- 
ter of the Milky Way. We adopt the lAU stan- 
dard values for the Galactocentric distance of the 
Sun, i?o = 8.5 kpc, and the circular rotational 
velocity of the LSR at the distance of the Sun, 
Oo = 220 km s^^. We also assume a flat Galactic 
rotation curve around the Sun. For almost all the 
parts, we adopt a Solar motion of 20 km s~^ rela- 
tive to the LSR to the direction of (aigoo, <5i9oo) = 
(18^^, 30°) for the radial velocity (radio definition 



of the LSR) . Hereafter this Solar motion is referred 
to as the standard Solar motion, and the LSR with 
the standard Solar motion is referred to as the tra- 
ditional LSR. In the following analyses, we adopt 
the systemic velocity of S Per of —38.5 km s~^ 
(Diamond et al. 1987) in the traditional LSR. On 
the other hand, because the Galactic parameters 
are known to involve some uncertainties (McMil- 
lan & Binncy 2010), we discuss the observational 
results adopting more realistic values for the So- 
lar motion in Section 5.2, which was determined 
from the average motion of a number of stars in 
the solar neighborhood observed with the Hippar- 
cos satellite. For the proper motion, we adopt the 
heliocentric coordinates. 

2. Observations 

Figure 1 shows sky positions of the observed 
sources. VLBI phase-referencing observations of 

5 Per at 22.2 GHz have been conducted together 
with a closely located continuum source as a posi- 
tional reference at a separation angle of 0.2°. This 
positional reference source was initially selected 
from the 6-cm northern sky catalogue (Becker et 
al. 1991) and the Texas survey (Douglas et al. 
1996). Later we found that this source is cata- 
logued under the number 143 in the 21-cm radio 
continuum survey of the Galactic plane by Kallas 

6 Reich (1980) (the so-called "KR survey") and 
an outer Galaxy VLA survey by Fich (1986) based 
on the KR survey, so that the source is referred to 
hereafter as KR143. Imai (2001) for the first time 
confirmed that KR143 is detectable with VLBI 
observations at 22 GHz, while we could not find 
an optical counterpart for KR143. We have con- 
ducted a total of eight-epoch VLBI observations 
over six years with the ten VLBA antennas. The 
observing epochs are listed in Table 1. 

Two intermediate-frequency (IF) bands with 
16-MHz bandwidth in left-hand circular polariza- 
tion were recorded. These IF bands are continuous 
in frequency. The H2O maser emissions of S Per 
has been received with the lower IF band (IFl) 
while the upper IF band (IF2) was used to improve 
the Signal-to-Noise Ratio (SNR) for KR143. We 
have also conducted phase-referencing observa- 
tions of a pair of KR143 and ICRF 0244+624 with 
a separation angle of 4.8° in the last five epochs 
to check the position stability of KR143 with re- 



spect to the positionally reliable ICRF source. A 
larger separation angle will take a longer telescope 
slew time, and therefore a longer switching duty 
cycle. In general, the larger the separation angle 
and/or the longer the switching duty cycle is, the 
more difficult is phase-referencing at 22 GHz due 
to stochastic tropospheric phase fluctuations. We 
selected a lower observing frequency of 15.3 GHz 
for the pair of KR143 and ICRF 0244+624 to re- 
duce a coherence loss effect due to the tropospheric 
phase fluctuations. The switching cycle time, tg, 
for the phase-referencing is listed in Table 1. 

In the first three epochs, 20-min phase-referencing 
sequence at 22 GHz was interleaved by 2- to 3-min 
calibrator observations. In the last five epochs, 20- 
min phase-referencing sequences at 22 GHz and 
15 GHz were conducted in turn, and each se- 
quence was interleaved by short calibrator scans 
at the same frequency. VLBI cross-correlation 
was carried out at the VLBA correlation center in 
Socorro to generate VLBI fringe data with 1024- 
frequency channels, corresponding to a velocity 
resolution of 0.2107 km s"-'^. In the first-epoch 
cross correlation, a wrong position with an offset 
by 4 arcsec was used as a phase tracking center of 
S Per. This mistake was recovered in the course 
of data reduction process. The observed pairs of 
sources in our monitoring program are listed in 
Table 2 together with their separation angles. 

3. Data reduction 

We analyzed the VLBI data using the stan- 
dard NRAO data reduction software. Astronom- 
ical Image Processing Software (AIPS), version 
31DEC06. There were two data reduction series: 
one for the pair S Per and KR143, and another 
one for the pair KR143 and ICRF 0244+624. 

3.1. Initial calibrations 

Initial calibration tasks were common for both 
the pairs of sources. First, the amplitude calibra- 
tion was conducted by using the measured system 
noise temperatures and gain variation data. In the 
next step. Earth Orientation Parameters (EOP) 
errors, two-bit sampling bias in the analogue-to- 
digital (A/D) conversion in VLBI signal process- 
ing, antenna parallactic angles, and ionospheric 
dispersive delays were corrected. 

In the first three epochs, NRAO 150, J0102+5824 



and 3C 84 have been observed as calibrators at 

22.2 GHz, while the former two have been ob- 
served as calibrators in the last five epochs. At 

15.3 GHz, NRAO 150 and J0102+5824 have been 
observed as calibrators in the last five epochs. 
Among them, NRAO 150 was used both in the 
data reduction series for calibrating relative clock 
offsets and phase offsets between the independent 
antennas. Phase differences between the two IF 
bands were also removed by using the NRAO 150 
fringe data. In the next steps of data analysis, a 
customized data reduction path was used for each 
pair of the sources. 

3.2. Data reduction for S Per and KR143 
at 22 GHz 

For the maser data analysis, the Doppler shift 
in S Per spectra due to the Earth rotation and 
Earth's Kepler motion were corrected with the 
AIPS spectrum channel correction task, CVEL. 
The observed frequencies of the maser lines were 
converted to radial velocities with respect to 
the traditional LSR using the rest frequency of 
22.235080 GHz for the H2O 616 - 523 transition. 
The spectrum deformation due to the receiver 
complex gain characteristics was calibrated with 
the NRAO 150 spectrum in the CVEL process. 

The 4 arcsec error in the position of S Per's 
phase tracking center in the first epoch was cor- 
rected with the AIPS visibility phase correction 
task UVFIX. Note that, although the S Per vis- 
ibilities might be subject to the coherence loss 
as a result of the cross correlation task with the 
phase tracking center error, the following astro- 
metric analysis is not affected. 

We inspected the cross-power spectra for all 
the baselines in order to select frequency chan- 
nels of emission which would be unresolved and 
not have rapid amplitude-time variation. We then 
carried out fringe-fitting by using the AIPS task 
FRING, with an option of no delay search for 
the selected channel. The fringe-fitted data were 
image-synthesized with a 50-/ias pixel for a region 
of 4096 X 4096 pixels, roughly a 200 x 200-square 
mas region by using the AIPS task IMAGR. Using 
this image as an initial model, phase and ampli- 
tude self-calibration was conducted. 

With the above processes from the fringe-fitting 
to the self-calibration, we obtained the complex 



gain calibration data, which were applied to the re- 
maining frequency channels to make image cubes 
for the same maser region. The images were pro- 
duced for all channels separately in the velocity 
range of —23 and —61 km s~^ with the step of 
0.2107 km s^^ in the traditional LSR. To pick up 
maser emission initially, we used the 2-D Gaussian 
component survey task SAD in each velocity chan- 
nel map. Based on the results of SAD, we visually 
inspected the surveyed Gaussian components one 
by one. We stored the image pixels larger than 3-(t 
noise level for a visually confirmed maser emission. 
Hereafter we define a maser "spot" as an emission 
in a single velocity channel, and a maser "feature" 
as a group of spots observed in at least two con- 
tinuous velocity channels at a coincident or very 
closely located positions. 

For KR143 imaging, the Doppler correction and 
bandpass calibration were skipped because the 
spectrum was averaged for the continuum source. 
Since there was no S Per's maser emission in IF2, 
the complex gain calibration data for IFl was ap- 
plied to IF2, whose phase and delay biases from 
IFl was already removed in the initial calibration. 
The synthesized image of KR143 with the pixel 
size of 50 /ias was obtained by using two IF bands, 
so that the image SNR is higher than that with 
the single IF band by a factor of \/2. The SNR 
of KR143 phase-referenced images ranges between 
18 and 37 at 22 GHz. 

3.3. Data reduction for KR143 and ICRF 0244+624 
at 15 GHz 

The Doppler and bandpass correction were 
skipped for both sources in processing of the 
15 GHz data as unnecessary ICRF 0244-^624 
was imaged with the 50-/xas pixel size to obtain 
the complex gain calibration data. This calibra- 
tion data were applied to the KR143 fringe at 
15 GHz before imaging. Synthesized images of 
KR143 were made with the 50-/zas pixel size. The 
SNR of KR143 phase-referenced images ranges 
between 7 and 17 at 15 GHz. 

4. Results 

4.1. KR143's position 

Phase-referenced radio maps of KR143 are 
dominated by a weak single component both at 15 
and 22 GHz. For the present astrometric study. 



we firstly determined the position of KR143 with 
respect to ICRF 0244+628 at 15 GHz, which had 
not been determined to high precision previously. 
In this paper, the position of ICRF 0244+624 
in J2000 is adopted as 2*^44^57". 696746 and 
+62°28'6".515030 in right ascension and decli- 
nation, respectively. Since ICRF 0244+624 has 
shown a very simple structure in the radio maps, 
we treated the peak position of the image as a po- 
sitionally stable point in the sky. We assume that 
KR143's position at 22 GHz is consistent with 
that at 15 GHz, and astrometric analysis of S Per 
was determined with respect to this position. 

Since images of KR143 have rather low SNR, a 
Gaussian fitting for the detected component may 
be affected by a rather severe image distortion due 
to a high image noise even if the object consists of 
only a single component. In addition, the phase- 
referencing image may be distorted by the resid- 
ual phase errors, introduced, e.g., by the static 
tropospheric delays. We used the peak position 
of KR143's image as the positional reference in 
order to minimize those effects in determining the 
position. The position measurement error was cal- 
culated from the beamwidth over the SNR at the 
peak. 

For the astrometric measurements of Galac- 
tic sources, a positional reference should be a 
background extragalactic source such as an Ac- 
tive Galactic Nuclei (AGN). However, because 
we hardly know whether KR143 is a suitable 
background source, we then investigated how 
positionally stable KR143 is with respect to 
ICRF 0244+624. Figure 2 shows a time variation 
of the KR143's position in the last five epochs 
of our experiment. Assuming that KR143 is a 
Galactic source for which the annual parallax is 
significant, we carried out a Levenberg-Marquardt 
non-linear least-square analysis to obtain the an- 
nual parallax and proper motion, but could not 
find a reliable solution because the y^ value per 
degree of freedom was quite large, no less than 57. 
This suggests that the annual parallax and proper 
motion are below our threshold because KR143 is 
sufficiently far from the Sun. This threshold is de- 
termined by the accuracy of the relative position 
measurements in phase-referencing for the pair of 
the sources. Figure 2 also displays a sinusoidal 
pattern of our least-square analysis result as well 
as another sinusoidal pattern for a source at the 



distance of 4 kpc, for comparison. 

The spatial distribution of the image peaks of 
KR143 at 15 GHz are shown in the top plot of 
Figure 3. This plot also shows the l-cr ellipse (a 
solid line). One possible reason for the observed 
change position of the peaks is inaccuracies in the 
VLBI correlator model for the phase-referencing 
astrometry for the pair of the observed sources. 
We conducted Monte Carlo phase-referencing ob- 
servation simulations for the pair of KR143 and 
ICRF 0244+624 to examine whether the observed 
position change could be explained by the un- 
certainties of the VLBI correlator model. For 
generating simulated phase-referencing fringes, we 
used a VLBI observation simulator, ARIS (Astro- 
nomical Radio Interferometer Simulator, Asaki et 
al. 2007b). In the simulations, we assume fiux 
densities of KR143 and ICRF 0244+624 of 20 
and 500 mJy, respectively, which were determined 
from our observations. The sources were assumed 
to be point-like sources. The simulated observa- 
tion schedule and observing system settings such 
as the recording bandwidth and A/D quantiza- 
tion level were adopted from actual VLBA obser- 
vations. There were two simulation series: one in- 
cluded l-cr static tropospheric excess path length 
error to the zenith (zenith EPL error) of 3 cm 
(e.g., Reid et al. 1999), and the other one used 
the zenith EPL error of 6 cm. We made 100 tri- 
als for each of the zenith EPL error cases. The 
simulation results are shown in the middle and 
bottom plots of Figure 3 for the l-cr zenith EPL 
error of 3 and 6 cm, respectively. In the case of a 
zenith EPL error of 6 cm, the simulation results 
arc consistent with the observation results. Since 
some of our observations were conducted during 
the summer period, it is quite natural that the 
zenith EPL error would be around 6 cm. We con- 
clude that the KR143's position change is mainly 
due to the zenith EPL error of several centimeters 
in the VLBI correlator model. Hereafter we refer 
a relative position error for a pair of point sources 
with the phase-referencing due to uncertainties in 
a VLBI correlator model for the atmospheric EPL, 
antenna positions, EOP, and other phase error fac- 
tors as an astrometric error. 

The measured flux densities of KR143 are 
shown in Figure 4 together with the data from 
the previous work by Douglas et al. (1996), Kallas 
& Reich (1980), Condon et al. (1998), Becker et 



al. (1991), Fich (1986), and Imai (2001). Note 
that, for the measurements of the flux density of 
KR143 from our radio images, we conducted a 
phase-self-cahbration process. Figure 4 resembles 
a synchrotron power-law spectrum, making it very 
likely that KR143 is an AGN. Ahhough there is 
currently no concrete evidence that KR143 is not 
located in or near our Galaxy because of the un- 
known redshift, we assume throughout this paper 
that KR143 is an AGN. The relative position of 
KR143 on the sky with respect to ICRF 0244-^624 
is determined with the standard errors of 61 and 
116 /xas in right ascension and declination, re- 
spectively, from all the five epochs. The absolute 
position of KR143 is listed in Table 3. 

4.2. Astrometry of S Per 

In the phase-referencing astrometric analysis 
for S Per, we investigated the time variation of 
the peak positions of the maser spots with re- 
spect to the peak position of the 22-GHz phase- 
referenced radio image of KR143. To estimate the 
astrometric error with the VLBI phase-referencing 
for the pair of S Per and KR143, we conducted 
Monte Carlo simulations of phase- referencing ob- 
servations for the pair of these sources at 22.2 GHz 
with ARIS. In the simulations, we assumed that 
flux densities of S Per and KR143 of 1 Jy and 
20 mjy, respectively, and that S Per was a contin- 
uum source. The tropospheric zenith EPL error of 
6 cm was set as found in the simulation described 
above for the pair of KR143 and ICRF 02444-624. 
Due to the small separation angle and fast antenna 
switching, the simulations show that the l-a errors 
in the relative position are 6 ^as and 14 /las in 
right ascension and declination, respectively. We 
used these standard deviations to provide the as- 
trometric error for the pair of the sources in the 
following analysis. 

Figure 5 shows the spatial distribution of the 
H2O maser emission of S Per for all the epochs. 
The maser emission is spread within a roughly 
circular structure with a radius of about 60 mas. 
We identified 77 maser features consisting of 478 
maser spots in all the available epochs. For ob- 
taining the trigonometric parallax and proper mo- 
tion of the star, we conducted the following four 
steps. Firstly, we selected a limited number of 
maser spots to evaluate the annual parallax. Our 
selection criteria arc as follows: (1) maser spots 



at epochs D, E, F, and G were all detected at the 
same velocity channel with the SNR of at least 
6; (2) in addition to the above epochs, they were 
detected with the SNR of at least 6 at least at 
one of the other epochs at the same velocity chan- 
nel; (3) no significant velocity drift larger than 
0.25 km s^^ yr^^ was observed for the features 
involving the selected spots; and (4) morphology 
does not significantly change between the epochs. 
We selected 44 maser spots from eight maser fea- 
tures. 

Secondly we carried out the Levcnberg-Marquardt 
least-square analysis for the i-th maser spot to ob- 
tain the initial position, (Aa^, A(5J^), at epoch A, 
the constant proper motion, (/i^cos(5, /i^), along 
with the annual parallax, tt', both in right ascen- 
sion and declination. In the initial trial, reduced 
X^ values of the analysis were in the range from 
a few tens to a few hundreds. These large values 
indicate that there must be unaccounted errors in 
the analysis. As explained above, we expect that 
the astrometric error in our observation is 6 and 
14 /xas for right ascension and declination, respec- 
tively. On the other hand, internal morphological 
time variations of the maser features during the 
monitoring period of over six years cannot be neg- 
ligible, as already pointed out in previous studies 
(e.g., Imai et al. 2007). For diagnostic purposes, 
we added position uncertainties of 0.2 mas for 
epochs A and H, 0.1 mas for epochs B and C, 
and 0.05 mas for epochs D, E, F, and G to the 
astrometric error both for right ascension and dec- 
lination. This resulted in making the reduced x^ 
values equal to unity in our least-square analy- 
sis. This additional position uncertainty which 
we attribute to the unpredicted time variations 
in maser morphology is referred to hereafter as 
a maser morphology uncertainty. Table 4 lists 
the least-square analysis results for each of the 
selected maser spots. Figure 6 shows the analy- 
sis result of one of the maser spots (spot-ID 39 
in Table 4). Figure 7 shows the histogram of 
the evaluated annual parallaxes for the 44 maser 
spots. The weighted average of the evaluated an- 
nual parallaxes is 0.414 ± 0.006 mas. 

We thirdly conducted a combined fitting with 
our least-square analysis for the 44 maser spots to 
obtain the trigonometric distance to S Per. In this 
step, we assumed that all the maser spots have 
a common annual parallax, tt* , which is one of 



the fitting parameters. This combined fitting re- 
sult for TT* is 0.413 mas, almost identical to the 
weighted average of the 44 annual parallaxes. The 
parameter fitting error in tt* was 0.006 mas from 
a covariance matrix of our least-square analysis. 

In obtaining the annual parallax, two system- 
atic errors have to be considered. One is the as- 
tromctric error, and the other one is the maser 
morphology uncertainty which provides a com- 
mon position offset to the member spots in a cer- 
tain maser feature. The most conservative esti- 
mate of the obtained annual parallax error be- 
comes 0.040 mas, which is calculated from the 
above statistical error of 0.006 mas multiplying 
by a factor of yll under an assumption that the 
position error of the selected maser spots are not 
independent. Hachisuka et al. (2009) applied a 
method to average positions of a number of maser 
spots from three features of W3(0H) water masers 
before their parallax fitting in order to avoid the 
maser morphology uncertainty, and achieved a 6- 
/^as accuracy in the annual parallax measurement. 
We note that the impact of the astrometric er- 
ror on the annual parallax measurement is not re- 
duced even if a number of maser features are simul- 
taneously processed. On the other hand, the im- 
pact of the maser morphology uncertainty can be 
reduced in processing multiple maser features be- 
cause the position offset due to this error becomes 
randomized. Since there are eight maser features 
in our analysis, the error quantity of 0.040 mas is 
over estimated. Figure 8 shows the position resid- 
uals of the selected 44 maser spots after removing 
the obtained annual parallax of 0.413 mas, proper 
motions and initial positions. It is noticed that, al- 
though a small positive trend can be seen in right 
ascension at epoch D, there are not apparently 
obvious systematic errors at central four epochs 
from D to G. It is also recognizable that system- 
atic errors are patially seen for specific maser fea- 
tures possibly due to the time variations in the 
maser morphology (see, for example, spot-ID 41- 
44 (feature-ID 8) at epochs D, E, and F). 

To evaluate contributions of those systematic 
errors to the obtained annual parallax, we con- 
ducted simulations of the combined fitting analy- 
sis by using fake data sets imitating the selected 
maser spots with the position measurement error, 
astrometric error, and maser morphology uncer- 
tainty. A random position offset due to the po- 



sition measurement error to each maser spot at 
each epoch. A position offset due to the astromet- 
ric error is applied to all the maser spots at each 
epoch. A position offset due to the maser mor- 
phology uncertainty is applied to all the maser 
spots included in each of the eight maser fea- 
tures at each epoch. We investigated the follow- 
ing error cases: (1) only position measurement 
error; (2) only astrometric error; (3) only maser 
morphology uncertainty; and (4) a combination 
of all the three errors above. We generated 100 
data sets for each of the four cases and carried 
out the Levenberg-Marquardt least-square analy- 
sis. Root-mean-square errors from the average of 
the 100 obtained annual parallaxes were 3, 9, 15, 
and 17 /xas for the cases of 1, 2, 3, and 4, respec- 
tively. The maser morphology uncertainty dom- 
inates the annual parallax error in our combined 
fitting while the astrometric error corresponds to a 
relatively small part. It is reasonable to adopt the 
value of 17 /xas as an upper limit of the trigono- 
metric distance uncertainty, corresponding to the 
trigonometric distance of 2.42lg;|!)g kpc. We then 
conducted our least-square analysis to obtain in- 
dividual proper motions of all the 478 maser spots 

with TT*. 

The obtained proper motions of the maser spots 
are considered to be a combination of an expand- 
ing flow of the circumstellar envelop (CSE) and the 
proper motion of the red supergiant in the Milky 
Way. We finally conducted a model fitting anal- 
ysis assuming a spherically expanding flow model 
described by Imai et al. (2000) and Imai et al. 
(2003). In the model fitting, we estimated the dis- 
tribution centroid indicating the stellar position as 
well as the stellar proper motion separated from 
the expanding shell motion of the CSE. We ob- 
tained the stellar proper motion in right ascension 
as /i* cos 5 = —0.49 ± 0.23 mas yr""'^, and in dec- 
lination as /i| = —1.19 ± 0.20 mas yr~^. This 
proper motion is consistent with that listed in the 
Tycho-2 Hipparcos catalogue, (—0.1 ± 1.8, — 2.8± 
1.8) mas yr^^, within its uncertainty (Hog et al. 
2000). The centroid of the maser distribution was 
determined with an accuracy of 8 mas in both right 
ascension and declination. Results of our astro- 
metric analysis for S Per are listed in Table 5. A 
residual time variation of the maser emission after 
removing the annual parallax and stellar proper 
motion is shown in Figure 9. Internal proper mo- 



tions of the identified 478 maser spots are shown 
in Figure 10 along with their radial velocities. 

5. Discussions 

5.1. Distance to S Per 

Our monitoring program determined the trigono- 
metric distance to S Per to be 2.42^QQg kpc with- 
out model assumptions. Marvel (1996) estimated 
a distance to this star as 2.3 ± 0.5 kpc by means 
of a statistical parallax method, in which relative 
proper motions of the H2O maser spots are fitted 
to an expanding shell where proper motions of the 
outermost spots are equal to the radial velocities 
of the maser spots with no significant proper mo- 
tions. The statistical parallax distance coincides 
with the trigonometric one within the measure- 
ment uncertainties. The trigonometric parallax 
greatly improved the accuracy of the distance: 
the uncertainty is five times better than that of 
the statistical parallax. 

S Per is considered to belong to Per OBI be- 
cause the radial velocity of S Per is consistent with 
that of this OB association (Bidelman 1947). The 
double cluster, h + x Per, whose distance was de- 
termined to be 2.34 kpc with the 2-percent accu- 
racy from the H-R diagram (Slesnick et al. 2002) 
is also included in Per OBI, so that S Per is in the 
vicinity of the double cluster in the OB associa- 
tion. The apparent radius of Per OBI is about 4°, 
corresponding to a ^^ 170-pc radius at the distance 
oi h + X Per. S Per is located on the far side of 
Per OBI. 

Recent VLBI astrometric observations of Galac- 
tic H2O and CH3OH maser sources in the Perseus 
arm have been used to measure their trigonomet- 
ric distances. Figure 11 shows the locations of 
such Galactic masers close to S Per in the Perseus 
arm in a face-on view of the Milky Way (Xu et 
al. 2006; Sato et al. 2008; Moellenbrock et al. 
2009; Moscadelh et al. 2009; Rcid et al. 2009a). 
In the figure, we depict only the H2O masers for 
NGC 281 West; measurements using the methanol 
masers result in a smaller trigonometric distance 
than using the H2O masers (Rygl et al. 2010). 
From this figure, it is clear that S Per, as well 
as NGC 281 West, are located at the outer edge 
of the Perseus arm. Assuming a pitch angle of 
the spiral arm of 16.5° (Reid et al. 2009b), we 
can estimate the width of the Perseus arm from 



the difference in locations between the outer edge 
maser sources and those located in the inner part 
of the arm. Hence, the estimated width of the 
Perseus arm is at least ^ 0.7 kpc. 

5.2. Three-dimensional motion of S Per 
in the Milky Way 

The stellar proper motion of S Per enables us 
to discuss the three-dimensional kinematics of the 
evolved star in the Milky Way. In this subsec- 
tion, we use realistic values of the Solar motion ^ , 
(Uq, Vq, Wq)= (7.5, 13.5, 6.8) km s"! (Francis & 
Anderson 2009; see also discussions in McMillan 
& Binney 2010). Let us estimate a non-circular 
motion of S Per in LSR at the red supergiant. 
This star has a rather large non-circular motion 
of 15 km s^^ at the Galactocentric distance of 
10.4 kpc: the non-circular velocity components 
are 3.3 ± 2.3 km s^^ toward the Galactic center, 
12.9 ± 2.9 km s^^ toward the anti-direction of the 
Galactic rotation, and 6.3±3.1 km s~^ toward the 
south pole. 

Recent VLBI astrometric observations of Galac- 
tic maser sources have revealed large non-circular 
motions of several objects in the Perseus arm. 
Figure 11 presents the non-circular velocity com- 
ponents of S Per along with those from other 
maser objects in the Perseus arm. Although this 
non-circular motion is not sensitive to Rq and Oq 
with the 10-percent variation, we have to admit 
that, if there is a 10-percent dip in the Galac- 
tic rotation curve at the Galactocentric distance 
of 9-10 kpc (Sofue et al. 2009), the lag of S Per 
from the Galactic rotation {V component) be- 
comes smaller. However, such a rotation curve 
cannot remove all the observed non-circular mo- 
tions of the maser sources in the Perseus arm, or 
even makes several of them larger. Moreover, such 
peculiar motions can be found for Galactic masers 
in star forming regions (SFRs) not only along the 
Perseus arm, but also along the local spur. Outer 
and Inner arms (Reid et al. 2009b; Rygl et al. 
2010). These two papers also discuss a large lag 
of the V component of the SFRs from a flat rota- 
tion curve together with a Hipparcos Solar motion 



^ Uq is the velocity component directed to the Galactic cen- 
ter, V0 is the component toward the circular rotational 
direction in the Galactic plane, and Wq is the component 
perpendicular to the Galactic plane and directed to the 
Galactic north pole. 



of Dehnen & Binney (1998). On the other hand, 
McMihan & Binney (2010) reanalyzed the sam- 
ples pubhshed by Reid et al. (2009b) in order to 
estimate the Galactic rotation as well as the Solar 
motion. Their fitting result of the Vq component 
of 11 km s^^, instead of 5.2 km s^^ of Dehnen 
& Binney (1998), can reduce such the large lags. 
However the Galactic SFRs still seems to have 
10 — 20 km s~^ peculiar motions even with the 
revised Solar motion, and those non-circular ve- 
locity components in the Galactic plane somehow 
show randomness. These non-circular motions do 
not quantitatively match a standard understand- 
ing of the stationary density-wave theory (Lin & 
Shu 1964) and a standing galactic shock solution 
in a tight-winding spiral potential (Roberts 1969; 
Roberts & Yuan 1970). 

One of the reasons for such peculiar motions 
may be attributed to an intrinsic space motion of 
each SFR. Lee & Lim (2008) found that there are 
new SFRs in the southern Hi shell in Per OBI, 
and that the new stars have larger southern veloc- 
ities due to the up-to-date cloud shuffling caused 
by ^ + X Per. According to the cloud-shuffling 
model, a molecular cloud is swept up by UV ra- 
diation pressure from nearby OB stars, perhaps 
together with supernova explosions. Star forma- 
tions are triggered in the cloud to form OB asso- 
ciations. At the same time, the compressed cloud 
will be accelerated in the anti-direction to the orig- 
inal sources causing the cloud shuffling. Cappa & 
Herbstmeier (2000) reported on an Hi superbub- 
ble with 2°-radius centered a.t h + x Per. How- 
ever, S Per, located inside the superbubble, does 
not move away from the double cluster, so that 
this star does not seem to be formed in a process 
of the cloud shuffling caused hy h + x Per. 

Although we cannot find a trace of cluster shuf- 
fling sources other than h + x Per at this moment, 
the cloud shuffling model still provides us an at- 
tractive idea in order to explain the large non- 
circular motion of the OB association and S Per. 
It is interesting that Lee & Lim (2008) also re- 
ported that proper motions of stars in Per OBI 
based on the Hipparcos measurements show a ve- 
locity component downward to the Galactic plane 
(toward the South Galactic pole), so that the non- 
circular motion of S Per may be representative of 
a motion of the molecular cloud in which the OB 
stars and S Per were formed. The plausible sce- 



nario is as follows: Considering the age oi h+x Per 
is 12.8 ± 1.0 Myrs (Slesnick et al. 2002), a cloud 
shuffling took place near the precursor molecular 
cloud of Per OBI around a few tens of Myrs ago 
because of radiation pressure from precedent OB 
stars and/or supernova explosions in the Galactic 
plane. The molecular cloud later forming Per OBI 
has been accelerated in the out-of-disk direction. 
S Per was formed in the molecular cloud ahead of 
a large SFR around the double cluster. This cloud 
shuffling also might affect the in-disk non-circular 
motions of the cloud. Another item of observa- 
tional evidence for a cloud shuffling causing the 
non-circular motion was provided by Sato et al. 
(2008) and Rygl et al. (2010) for NGC 281 West. 
It was suggested that large non-circular motions of 
NGC 281 West masers were due to the supershell 
generated by supernova explosions. 

Another new idea to explain the non-circular 
motions of the Galactic SFRs arises from large 
N-body/hydrodynamics numerical simulations. 
Baba et al. (2009) showed that, with their high- 
resolution numerical model of a spiral galaxy 
and their original N-body/gas simulation code 
(Saitoh et al. 2009), the non-circular motions 
of the SFRs were explained by natural conse- 
quence of non-linear interactions between non- 
stationary spiral arms and multi-phase interstel- 
lar media. With their numerical simulations of 
a disk galaxy, they found that young stars with 
ages < 50 Myr in a spiral galaxy have large and 
random non-circular motions. This idea can be 
applied to S Per and Per OBI because the ages 
of those objects are younger than 50 Myr. The 
N-body-f-hydrodynamics simulations also revealed 
that molecular clouds in a spiral galaxy can have 
a velocity dispersion of several km s~^ in the ver- 
tical direction to the disk (J. Baba, 2009, private 
communications), so that a rather large peculiar 
motion of Per OBI vertical to the Galactic plane 
could be explained without a cloud shuffling. Since 
we have quite limited samples with the measur- 
able non-circular motions of the young objects in 
the Milky Way, it is very important to measure 
the distances and proper motions of the spiral 
arm objects to compare with those simulation re- 
sults. A deep understanding of the dynamics of 
the Galactic spiral arms and the origin of the non- 
circular motions of SFRs will be possible in near 
future when massive numerical simulations can 



be combined with a much larger amount of VLBI 
astrometry samples. 



6. 



Conclusions 



We have conducted phase-referencing VLBI 
monitoring for the Galactic red supergiant, S Per. 
The distance to S Per was determined to be 



2.42 



+0.11 
-0.09 



kpc based on the trigonometric paral- 
lax. Together with the stellar radial velocity and 
obtained proper motion, we investigated the three- 
dimensional motion of S Per in the Milky Way, 
showing the non-circular motion of 15 km s~^ from 
the Galactic circular rotation, which is mainly 
dominated by the anti rotation direction com- 
ponent of 12.9 ± 2.9 km s^^. This non-circular 
motion may be representative of the motion of the 
young star cluster. Per OBI, in the Milky Way. 
Recent large N-body/hydrodynamics numerical 
simulations show that the non-circular motions of 
star forming regions are not exceptional but can be 
seen everywhere in a spiral galaxy. The cloud shuf- 
fling is another plausible mechanism to give such a 
large non-circular motion to the Per OBI associa- 
tion. To reveal these non-circular motions as well 
as the other maser sources on the spiral arms, a lot 
of sample points are needed. In the near future, we 
will be able to obtain a set of trigonometric mea- 
surements with the three-dimensional motions of 
the maser sources with phase-referencing VLBI. 
We will be able to investigate the galactic struc- 
ture along with its kinematics with a support of 
the large-scale N-body/hydrodynamics numerical 
simulations. 
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proper motions, and initial positions. The abscissa represents the maser spot ID as listed in Table 4, and 
the ordinate is the position residuals in mas. For each epoch, the position residuals in right ascension 
and declination are shown in the upper and lower part, respectively. Error bars represent the position 
measurement error. 
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in the figure (1 to 8) indicate the maser feature-ID which was used for obtaining the stellar annual parallax. 
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Fig. 11. — Spatial positions and non-circular velocity components of S Per and other maser sources in 
the Perseus arm projected on the Galactic plane. The depicted sources are W3 (OH) (Xu et al. 2006), 
NGC 281 West (Sato et al. 2008), IRAS 00420+5530 (Moellenbrock et al. 2009), S 252 (Reid et al. 2009a), 
NGC 7538 (Moscadelli et al. 2009), and S Per (present paper). A dark-gray dashed curve indicates the Perseus 
arm obtained from Hh survey (Georgelin & Georgelin 1976) while a light-gray dashed curve indicates the one 
obtained from CO survey (Nakanishi & Sofue 2006). Those traces are adjusted to those for Rq — 8.5 kpc. 
Dotted curves indicate equi-distances from the Galactic center (8, 10, and 12 kpc). 
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Table 1: Observation epochs and duty-cycle times for VLBI phase- referencing. 









ts at 22 GHz 


ts at 15 GHz 


Epoch 


Date 


Time range 


(KR143 - 


(KR143 - 






(UTC) 


SPer) 


ICRF 0244+624) 


A 


2000 Nov 21 


03:10-08:10 


40 s 




B 


2002 May 3 


16:26-21:26 


40 s 


- 


C 


2002 May 29 


14:44 - 19:44 


40 s 


- 


D 


2003 Jul 13 


10:48 - 17:48 


35 s 


60s 


E 


2003 Oct 7 


05:10 - 12:09 


35 s 


60s 


F 


2004 May 2 


15:29-22:27 


35 s 


60s 


G 


2004 Aug 11 


08:46 - 15:46 


35 s 


60s 


H 


2007 Feb 22 


20:03 - 03:01 


35 s 


60s 
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Table 2: Observed pairs of sources in our monitoring program. 



Target Phase-referencing Separation Frequency Epoch 

cahbrator angle 

KR143^ SPct 0^2^ 22.2 GHz A-H 

KR143 ICRF 0244+624^ 4.8° 15.3 GHz D - H 



''KR143 is a positional reference of S Per in VLBI phase-referencing at 22 GHz. 

''ICRF 0244-1-624 is a positional reference of KR143 in VLBI phase-referencing at 15 GHz. 
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Table 3: Determined KR143's positions. 



Source RA (J2000) Dec (J2000) 

KR143^ 2'^22'"33". 520836 ± O^OOOOIS +58°48'13". 944455 ± 0".000256 



^Thc positional reference in the International Celestial Reference Frame (ICRF) is ICRF 0244+624. The position error of 
ICRF 0244+624 in ICRF (0.124 and 0.228 mas in right ascension and declination, respectively) is taken into consideration. 
This position is located at the origin of Figure 2 and of the top plot in Figure 3. 
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Table 4 
Proper motions and parallaxes for the S Per H2O maser spots. 



Spot ID 


AttA 


ASa 


Vlsr 


TT 


lla COS 5 


^J■5 


(Feature 
ID) 


[mas] 


[mas] 


[km s-i] 


[mas] 


[mas yr^"'"] 


[mas yr~^] 


1(1) 


-27.85 ±0.11 


0.11 ±0.03 


-33.58 


0.417 ±0.047 


-0.95 ±0.04 


-1.40 ±0.01 


2(1) 


-27.86 ±0.10 


0.26 ±0.03 


-33.79 


0.398 ±0.041 


-0.93 ±0.03 


-1.43 ±0.01 


3(1) 


-27.86 ±0.10 


0.28 ±0.03 


-34.00 


0.398 ±0.040 


-0.93 ±0.03 


-1.44 ±0.01 


4(1) 


-27.89 ±0.10 


0.27 ±0.03 


-34.21 


0.425 ±0.042 


-0.92 ±0.03 


-1.42 ±0.01 


5(1) 


-27.98 ±0.12 


0.48 ±0.04 


-34.42 


0.410 ±0.051 


-0.89 ±0.04 


-1.50 ±0.01 


6(2) 


-17.31 ±0.17 


-7.16 ±0.03 


-45.59 


0.403 ±0.046 


-1.35 ±0.05 


-1.43 ±0.01 


7(3) 


-46.82 ±0.10 


-27.82 ±0.03 


-29.57 


0.412 ±0.039 


-1.23 ±0.03 


-1.66 ±0.01 


8(3) 


-45.69 ±0.12 


-26.81 ±0.03 


-29.57 


0.454 ±0.039 


-1.37 ±0.04 


-1.73 ±0.01 


9(3) 


-46.83 ±0.10 


-27.79 ±0.03 


-29.78 


0.429 ±0.040 


-1.21 ±0.03 


-1.66 ±0.01 


10(3) 


-45.84 ±0.10 


-26.86 ±0.03 


-29.78 


0.470 ±0.038 


-1.32 ±0.03 


-1.71 ±0.01 


11(3) 


-45.46 ±0.09 


-26.52 ±0.03 


-30.21 


0.352 ±0.043 


-1.37±0.03 


-1.79 ±0.01 


12(4) 


-14.57 ±0.15 


-7.89 ±0.04 


-47.06 


0.475 ±0.052 


-0.81 ±0.05 


-1.76 ±0.01 


13(4) 


-13.52 ±0.14 


-8.53 ±0.03 


-47.69 


0.431 ±0.041 


-1.21 ±0.04 


-1.50 ±0.01 


14(4) 


-13.24 ±0.14 


-8.37 ±0.03 


-48.12 


0.411 ±0.042 


-1.30 ±0.04 


-1.54 ±0.01 


15(4) 


-12.08 ±0.11 


-8.78 ±0.03 


-48.12 


0.372 ±0.039 


-0.88 ±0.04 


-1.72 ±0.01 


16(4) 


-13.28 ±0.14 


-8.37±0.03 


-48.33 


0.412 ±0.043 


-1.29 ±0.04 


-1.54 ±0.01 


17(4) 


-11.36 ±0.13 


-10.00 ±0.03 


-48.33 


0.449 ± 0.042 


-0.90 ±0.04 


-1.54±0.01 


18(4) 


-7.85 ±0.12 


-9.12 ±0.03 


-49.38 


0.404 ±0.043 


-1.26 ±0.04 


-1.72 ±0.01 


19(4) 


-7.81 ±0.11 


-8.95 ±0.03 


-49.59 


0.386 ±0.043 


-1.27 ±0.04 


-1.77 ±0.01 


20(5) 


41.93 ±0.10 


-26.14 ±0.04 


-39.90 


0.384 ±0.051 


0.28 ±0.03 


-1.57±0.01 


21(5) 


41.93 ±0.10 


-26.15 ±0.03 


-40.11 


0.403 ± 0.043 


0.27 ±0.03 


-1.56 ±0.01 


22(5) 


41.88 ±0.10 


-26.15 ±0.03 


-40.32 


0.425 ±0.042 


0.27 ±0.03 


-1.54 ±0.01 


23(5) 


41.83 ±0.10 


-26.15 ±0.03 


-40.53 


0.401 ±0.045 


0.28 ±0.03 


-1.56 ±0.01 


24(6) 


20.53 ±0.10 


-43.94 ±0.03 


-43.48 


0.472 ±0.047 


-0.15 ±0.03 


-1.91 ±0.01 


25(6) 


20.29 ±0.10 


-44.08 ±0.03 


-43.69 


0.332 ±0.042 


-0.08 ±0.03 


-1.91 ±0.01 


26(6) 


20.04 ±0.10 


-44.47 ±0.03 


-43.90 


0.405 ±0.041 


-0.03 ±0.03 


-1.82 ±0.01 


27(6) 


20.12 ±0.12 


-44.67 ±0.03 


-44.11 


0.425 ±0.041 


-0.07 ±0.04 


-1.77±0.01 


28(6) 


20.28 ±0.14 


-44.78 ±0.03 


-44.32 


0.419 ±0.042 


-0.13 ±0.04 


-1.75 ±0.01 


29(6) 


20.27 ±0.12 


-44.75 ±0.03 


-44.53 


0.426 ±0.041 


-0.13 ±0.04 


-1.78 ±0.01 


30(6) 


19.21 ±0.10 


-44.09 ±0.03 


-44.74 


0.425 ±0.040 


-0.05 ±0.03 


-1.72 ±0.01 


31(6) 


20.26 ±0.13 


-44.81 ±0.03 


-44.74 


0.418 ±0.044 


-0.14 ±0.04 


-1.78 ±0.01 


32(6) 


19.10 ±0.10 


-43.90 ±0.03 


-44.96 


0.417 ±0.039 


-0.04 ±0.03 


-1.78 ±0.01 


33(6) 


19.05 ±0.10 


-43.97 ±0.03 


-45.17 


0.408 ±0.039 


-0.03 ±0.03 


-1.75 ±0.01 


34(6) 


19.03 ±0.10 


-44.06 ±0.03 


-45.38 


0.419 ±0.039 


-0.03 ±0.03 


-1.73 ±0.01 


35(6) 


19.03 ±0.10 


-44.14 ±0.03 


-45.59 


0.430 ±0.039 


-0.04 ±0.03 


-1.72 ±0.01 


36(6) 


19.06 ±0.10 


-44.22 ±0.03 


-45.80 


0.433 ±0.041 


-0.06 ±0.03 


-1.73 ±0.01 


37(6) 


19.11 ±0.10 


-44.32 ±0.04 


-46.01 


0.454 ±0.050 


-0.09 ±0.03 


-1.74 ±0.01 


38(7) 


28.43 ±0.09 


-49.64 ±0.03 


-41.79 


0.430 ±0.037 


0.02 ±0.03 


-1.93 ±0.01 


39(7) 


28.47 ±0.09 


-49.79 ±0.03 


-42.01 


0.430 ±0.037 


-0.00 ±0.03 


-1.89 ±0.01 
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Table 4 — Continued 






Spot ID 

(Feature 

ID) 


[mas] 


Ma 
[mas] 


Vlsr 
[km s-i] 


[mas] 


/Xc COS 5 

[mas yr~^] 


M.5 

[mas yr~-^] 


40(7) 
41(7) 
42(8) 
43(8) 
44(8) 


28.46 ±0.09 
28.48 ±0.09 
-9.84±0.15 
-9.79 ±0.11 
-9.80 ±0.11 


-49.80 ±0.03 

-49.83 ±0.03 

3.81 ±0.03 

3.91 ±0.03 

3.98 ±0.03 


-42.22 
-42.43 
-47.91 
-48.12 
-48.33 


0.422 ±0.037 
0.428 ±0.037 
0.404 ±0.041 
0.378 ±0.039 
0.351 ±0.039 


-0.01 ±0.03 
-0.04 ±0.03 
-0.67 ±0.05 
-0.68 ±0.03 
-0.67±0.03 


-1.89 ±0.01 
-1.87±0.01 
-1.03 ±0.01 
-1.06 ±0.01 
-1.07±0.01 
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Table 5: Astromctry analysis results of S Per. 



Annual parallax 


TT* ^ 0.413 ± 0.017 mas (2.4211^;^^ kpc) 


Position (J2000)^ 
(2000 Nov 21) 


Right Ascension: Declination: 
2'^22'"51^7106 ± 0^.0010 +58°35'11".444 ± 0".008 


Stellar proper motion 


[i*^co^6= -0.49 ±0.23 masyr-^ /i| = -1.19 ±0.20 mas yr"^ 



^Thc position errors of ICRF 0244+624 and KR143 are taken into consideration as well as the model fitting error described in 
Section 4.2. This position is located at the origin of S Per image of Figure 5. 
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